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In  this  paper  we  present  a  closed-form  solution  to  the  problem  of  optimally  charging  a  Li-ion  battery.  A 
combination  of  three  cost  functions  is  considered  as  the  objective  function:  time-to-charge  CITC),  energy 
losses  (EL),  and  a  temperature  rise  index  (TRi).  First,  we  consider  the  cost  function  of  the  optimization 
problem  as  a  weighted  sum  of  TTC  and  EL.  We  show  that  the  optimal  charging  strategy  in  this  case  is  the 
well-known  Constant  Current-Constant  Voltage  (CC-CV)  policy  with  the  value  of  the  current  in  the  CC 
stage  being  a  function  of  the  ratio  of  weighting  on  TTC  and  EL  and  of  the  resistance  of  the  battery.  Then, 
we  extend  the  cost  function  to  a  weighted  sum  of  TTC,  EL  and  TRI  and  derive  an  analytical  solution  for  the 
problem.  It  is  shown  that  the  analytical  solution  can  be  approximated  by  a  CC-CV  with  the  value  of 
current  in  the  CC  stage  being  a  function  of  ratio  of  weighting  on  TTC  and  EL,  resistance  of  the  battery  and 
the  effective  thermal  resistance. 

©  2015  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Battery  charging  is  a  problem  of  significant  interest,  especially  as 
the  battery-dependent  smart  devices  proliferate.  The  literature 
abounds  with  different  strategies  for  charging  batteries.  Among  the 
traditional  methods  of  charging,  the  simplest  is  the  constant  trickle 
current  charge  strategy,  which,  due  to  its  low  charging  current, 
requires  a  long  charging  time  (around  10  h)  [9];  constant  current 
strategy  with  higher  rates  of  current  requires  shorter  charging  time. 
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The  most  widely-used  traditional  strategy  is  the  constant-current 
constant-voltage  (CC— CV)  [9]  strategy,  in  which  a  constant  cur¬ 
rent  is  applied  to  the  battery  until  the  terminal  voltage  reaches  a 
specified  value,  and  afterward  the  charging  current  decreases  by 
applying  a  constant  voltage  to  the  terminals  of  the  battery.  In  Refs. 
[21,22],  a  multi-step  constant-current  charging  is  devised  for 
shortening  the  charging  time  and  prolonging  the  cycle  life  of  the 
battery.  Using  orthogonal  arrays,  Taguchi-based  methods  for  bat¬ 
tery  charging  [12,23]  present  a  systematic  method  to  find  the 
optimal  solution  with  guidelines  for  choosing  the  design  parame¬ 
ters.  In  Ref.  [13],  a  boost  charging  strategy  is  proposed  by  applying 
very  high  currents  to  close-to-fully  discharged  batteries.  In  pulse¬ 
charging  methods  [6,7,17,20,27],  the  battery  is  exposed  to  very 
short  rest  or  even  deliberate  discharging  periods  during  the 
charging  process.  Soft-computing  approaches  are  also  used  in  the 
optimization  of  battery  charging  profile.  In  Ref.  [24],  the  charging 
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problem  is  viewed  as  an  optimization  problem  with  the  objective 
function  of  maximizing  the  charge  within  30  min  using  a  multi¬ 
stage  constant  current  charging  algorithm  whose  optimal  solution 
is  obtained  via  an  ant-colony  approach.  In  Ref.  [25],  a  universal 
voltage  protocol  is  proposed  to  improve  charging  efficiency  and 
cycle  life  by  applying  a  charging  profile  depending  on  the  state-of- 
health  (SOH)  of  the  battery,  using  SOH  estimation  approaches  [26] 
in  the  optimization  process.  Recently,  in  Ref.  [11],  battery  charging 
is  considered  as  an  optimization  problem  with  cost  function  of 
time-to-charge  and  energy  loss  (as  we  do  in  this  paper),  but  they 
have  not  solved  the  problem  analytically;  rather  they  have  pre¬ 
sented  a  numerical  solution  to  the  problem.  Other  approaches,  such 
as  genetic  algorithm  and  neural  network  based  strategies  [16],  data 
mining  [2,10],  and  Grey-predicted  charging  system  [8]  have  also 
been  used  for  charging  batteries. 

In  this  paper  we  look  at  the  charging  problem  from  a  fresh 
perspective  using  optimal  control  theory,  and  our  goal  is  to  find  the 
optimal  current  profile  that  minimizes  a  specific  cost  function.  In 
this  sense,  different  objectives  may  be  embedded  in  the  cost 
function.  One  obvious  cost  function  is  the  time-to-charge  (TTC).  We 
prefer  to  minimize  the  charging  time  as  much  as  possible,  as  TTC 
reduction  contributes  to  user  satisfaction.  Another  important 
objective  is  the  energy  loss  (EL)  during  charging.  Reducing  the 
energy  loss  increases  the  charging  efficiency.  In  this  paper,  first  we 
use  an  integrated  cost  function  that  includes  both  the  TTC  and  EL. 
Then,  we  also  include  the  effects  of  temperature  into  account,  and 
the  cost  function  is  selected  as  a  linear  combination  of  three 
criteria:  time-to-charge,  energy  loss,  and  temperature  rise  index 
(TRI).  In  both  cases,  analytical  solutions  of  the  optimal  charging 
problem  are  derived. 

The  paper  is  organized  as  follows.  In  Section  2,  we  derive  an 
analytical  solution  for  the  optimal  charging  current  profile  to 
minimize  TTC  and  EL.  In  Section  3,  we  extend  this  approach  to  the 
case  where  temperature  rise  is  considered  as  well.  Section  4  is 
devoted  to  simulation  results  and  finally  we  conclude  the  paper  in 
Section  5. 

2.  Analytic  solution  for  optimal  charging  current  profile 

We  consider  a  simplified  equivalent  circuit  model  of  the  battery 
as  shown  in  Fig.  1.  The  theory  extends  naturally  to  more  complex 
models  involving  parallel  RC  elements  (shown  in  Fig.  2),  but 
analytical  closed  form  solutions  are  not  possible  in  the  latter  case. 
The  model  consists  of  a  voltage  source  corresponding  to  the  open- 
circuit  voltage  (OCV),  which  is  dependent  on  the  state  of  charge 
(SOC),  and  a  resistance  R0.  The  SOC  and  OCV,  are  represented 
respectively  by  s  and  V0.  The  OCV  is  a  nonlinear  function  of  SOC  and 
is  denoted  by  V0(s[k]). 

The  state  of  charge  is  zero  when  the  battery  is  totally  discharged 
and  it  is  one  if  it  is  completely  charged.  The  sampling  time  is 


R0  i  [/c] 


Ri 


denoted  by  A  (in  seconds).  We  assume  that  the  initial  and  final  SOC 
are  known:  s[0]  =  s0,  s[fcy]  =  sk  ,  where  A  is  the  charging  time.  We 
also  assume  that  the  maximum  allowed  value  of  the  terminal 
charging  voltage  is  vc,  that  is,  ijk]  <  vc  for  all  k. 

The  SOC  dynamics  for  the  battery  considering  the  foregoing 
model  are  as  follows: 

s[k  +  1]  =  s[k]  +  chi[k]  (1) 

where  ch  (in  1  /Amperes)  is  the  parameter  in  Coulomb  counting, 
given  by 

Ch  =  3600Cbatt  (2) 

where  Cbatt  (in  Ah)  is  the  battery  capacity,  assumed  to  be  known. 

Let  the  objective  function  be  a  combination  of  TTC  and  EL.  In 
other  words, 

kf-i 

JtE  =  vvt/t  +  We/e  =  Wtkf\  +  wE^  Roi2[k}&  (3) 

fc=o 

where  Jt  is  the  TTC  cost  function,  JE  is  the  EL  cost  function;  vvt  and 
wE  are  weights  on  the  TTC  and  EL  cost  functions,  respectively.  The 
resistance  of  the  battery,  i.e.,  R0,  is  assumed  to  be  known. 

The  charging  problem  then  could  be  formulated  as  follows: 
Minimize  jtE  subject  to: 


s[k  +  1]  =  s[k]  +  chi[k\  s[0]  =  s0  s  [fy]  =  skf 

(4) 

V0  (s[k[)  +  R0i[k]  <  vmax 

(5) 

i[l<\  <  imax 

(6) 

It  is  important  to  note  that  only  the  ratio  of  weights  affects  the 
optimal  current  profile  of  i [/<■].  Therefore,  by  dividing  (3)  by  wE,  we 
redefine  the  cost  function  as  follows: 


/  kf 

JtE  =]tE /we  =  Pt/t  +Je  =  Pt^fA  +  ^  Roi2[k]A  (7) 

k=0 

where  pt  =  wt/wE.  Also  note  that  when  the  current  is  injected  into 
the  battery,  the  V0  starts  to  increase  and  this,  in  turn,  causes  the 
terminal  voltage  to  rise,  until  it  reaches  i'max.  which  is  the 
maximum  allowed  terminal  voltage.  During  the  whole  charging 
process  the  current  should  not  exceed  imax,  which  is  the  maximum 
allowed  charging  current.  In  this  paper,  we  use  vc  for  i/max,  where  vc 
is  the  voltage  corresponding  to  SOC  of  1 ;  that  is 

vc  =  V, o(l)  (8) 
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Assume  that  at  time  fc] ,  the  terminal  voltage  v[ki]  reaches  vc  and 
let  us  denote  the  state  of  charge  at  time  fc i  as  Si .  After  time  k] ,  the 
terminal  voltage  should  be  fixed  at  the  constant  voltage  (CV)  vc\ 
hence,  for  fc  =  k] ,  kj  +  1 , . . . , kf  —  1,  the  dynamics  of  the  system  are 
as  follows: 


i[k]=^-(vc-V0(s[k])) 

(9) 

s[/<+  1]  =  s[k]  +  chi[k] 

(10) 

s[ki]  =  s-,  s^kf^s^ 

(11) 

Before  going  further,  let  us  define  a  new  equivalent  problem  as 
follows: 

Minimize 


9H[k]  _ 

m 


m 


dH[k\ 
ds  [k] 


s[k  + 1]  = 


3H[k] 


6A[fc  +  1] 

From  (17)  we  have 

f  [k]  -  -Ch^1]  k  =  0, 1,  ...ki  —  1 
From  (18),  we  can  write 
A[k]=A[k  +  l]  k  =  ki  -  1 , . . . ,  0  A[k-|]  =  v 


(17) 

(18) 

(19) 

(20) 

(21) 


/<i-i 

J tE  =  PtJt  +Je  =  Pth&  +  ^2  i?oi2[k]A  (12) 

k= 0 

subject  to: 


where  v  is  the  Lagrange  multiplier  associated  with  the  constraint 
s[k1]=Si.  Equation  (21)  implies  that  all  co-states  are  equal; 
therefore,  we  can  write 

A[k]  =  v  k  =  0,1, _ ki  (22) 


s[k  +  1]  =  s[k]  +  chi[k ]  s[0]  =  s0  s[/<i]  =  S! 


(13)  Based  on  (22),  equation  (20)  can  be  written  as 


This  problem  is  in  fact  the  minimization  in  the  stage  where  the 
terminal  voltage  is  below  vc  and  therefore  here  the  condition 
V0  (s[fc[)  +  R(p[k]  <  vc  is  not  shown  as  we  know  that  it  holds. 

Inspired  by  Refs.  [18]  and  [1],  we  solve  the  problem  in  three 
steps  as  described  below: 

O  Given  k ]  (when  the  terminal  voltage  constraint  becomes  active), 
find  the  optimal  current  profile  that  minimizes  the  energy  los¬ 
ses,  and  calculate  the  corresponding  energy  losses  as  a  function 
of  fc] . 

©  Generate  a  new  equivalent  cost  function  JtE<  consisting  of  the 
weighted  TTC  plus  the  fc]  -dependent  minimum  energy  loss 
obtained  in  step  1,  and  find  the  optimal  fc ]  based  on  this  cost 
function. 

©  Given  the  optimal  k ]  from  step  2,  evaluate  the  optimal  current 
obtained  in  step  1. 

In  the  first  step,  assuming  Iq  is  known,  we  find  the  optimal 
current  i*  [fc|/<i  ]  that  minimizes  the  energy  loss.  Having  this  optimal 
current  profile,  we  can  calculate  the  minimum  EL  cost  function 
_/E(fc]),  which  is  a  function  of  fc] .  In  the  second  step,  we  use  the 
partially  optimized  cost  function  JtE*  =  ptfc1A  +  _/E(fc])  and  we  find 
the  optimum  value  for  kh  say  fc).  In  the  third  step,  we  insert  the 
optimal  final  time  fcj  into  the  current  i*[fc|fc1]  (obtained  in  step  1)  to 
find  the  optimal  current  i*[k].  The  final  result  will  be J/gjfcj). 

The  first  stage  is  formulated  as  follows: 

Minimize 


i*[k]  =  _2^ok  k  =  0-1-->ki-1  (23) 

Note  that  equation  (23)  states  that  the  optimal  current  is  con¬ 
stant.  From  (19),  we  can  write 

s[k  +  l]  =  s[k]  +  chi[k]  (24) 

which  is  actually  the  dynamics  of  the  system.  Knowing  the  initial 
state  of  charge  (s0),  and  noting  the  optimal  current  in  (23)  is  con¬ 
stant,  we  have 


k  1  hc2 

s[k]  =  s0  +  chJ2  i[l\  =  s0  - 

2R° A 

Since  for  fc  =  k i  we  have  s[ki]  =  Si,  therefore 
ki  c\v 

Sl=S°-2 rOk 

Solving  for  v,  we  have 

2i?0A(s1  -  s0) 
k>  ch 

Inserting  (27)  into  (23),  we  have 

=  k  =  0,l,...,k1-l 

Mch 


(25) 


(26) 


(27) 


(28) 


JE(ki)  =  E  ^0i2  [k]  A 
k= 0 

subject  to 

s[k  +  1]  =  s[k]  +  chi[k]  s[0]  =  s0  s[ki]  = 

The  Hamiltonian  function  for  this  problem  is 

H[k]  =  i?0i2[k]A  +  A[fc  +  l](s[k]  +  chi[li}) 


(14) 


Inserting  28)  into  the  cost  function,  the  optimal  cost  function, 
given  fc i  is: 


(15) 


*i-i 

jm)  =  J2  R° 

k= 0 


Sjj-Sq 

l<  \Ch 


A  = 


RqA(Si  -  Sp)^ 

fci  cl 


(29) 


Now,  consider  step  2  and  define  the  cost  function  as 


(16) 


JtE*  =  Pth  A  +Mki)  =  Pth  ^  + 


KoA(st  -  s0) 

hc2h 


(30) 


The  following  equations  must  hold  for  the  optimal  solution  [5]: 


To  find  the  optimum  fci ,  the  following  relations  should  hold: 
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JtE*  (^1  “  1)  —  JtE*  (^1 ) 

(31) 

JtE.(ki+l)>JtE.(ki) 

(32) 

Inserting  (30)  into  (31)  and  (32)  we  obtain  two  second-order 
equations  in  term  of  Jq .  Solving  these  equations,  we  get  kf  and 
k|,  respectively,  for  relations  (31)  and  (32). 


±1  + 


1  i  4Ro(Si—  So) 

1  H  T72. - 


PK 


(33) 


The  optimum  Iq  is  ceil(kj  )  or  floor(fcf ).  Since  /q  -  k|  =  1,  we 
have  ceil(k^ )  =  floor(k| )  =  round((fcj"  +  k|)/2).  Thus, 


k\  =  round 


1  Rq(si  -s0)2 
4 +  Ord 


(34) 


A  more  convenient  way  is  to  treat  fe1  in  (30)  as  a  continuous 
variable  and  take  derivative  of  (30)  with  respect  to  k ^  as  follows: 


9JtE*(^i)  _  *  ftoA(si“so)2 

Wi 


(35) 


k\ 


si  ~  s0 
Ch 


(36) 


Note  that  if  we  neglect  1  /4  in  (34),  the  argument  of  the  rounded 
function  in  (34)  is  exactly  the  same  as  the  one  in  (36). 

Step  3  involves  inserting  (36)  into  (28)  to  find  the  optimum 
current 


simplification  of  the  dynamics  of  the  system  into  the  condensed 
condition  of  (40)  will  be  useful  in  the  next  section  where  we  derive 
an  analytical  solution  when  the  cost  function  includes  the  sum¬ 
mation  of  temperature  rises  as  well. 

It  should  be  noted  that  the  practical  meaning  of  the  parameters 
of  optimization  problem  (e.g.,  wt  and  wE  in  (3)  and  pt  in  (12))  is  to 
use  them  in  an  iterative  design  procedure  to  reach  the  desired 
performance.  For  example,  if  the  maximum  allowed  energy  loss  is 
£max  and  the  maximum  acceptable  time-to-charge  is  TTCmax,  then 
in  the  design  procedure,  wt  and  wE  should  be  selected  inversely 
proportional  to  TTCm3X  and  £max,  respectively;  that  is  vvt  °c  1  /TTCmax, 
wE  °c  1  /£max  and  then  iterate.  Or  equivalently,  pt  should  be  selected 
proportional  to  Emax/TTCmax',  that  is  pt°cEmax/TTCmax  and  then 
iterate  on  the  proportionality  factor.  As  £max/E£Cmax  is  approxi¬ 
mately,  the  allowable  average  power  loss,  pt  should  be  selected 
proportional  to  allowable  average  power  loss. 

3.  Optimal  charging  problem  considering  temperature 

In  this  section,  we  will  extend  the  cost  function  to  include  the 
battery  temperature  via  temperature  rise  index  (TRI,  to  be  defined) 
as  well  as  TTC  and  EL.  To  this  end,  we  need  a  temperature  model  for 
the  battery.  Refs  15]  and  [19]  describe  the  temperature  model  of 
the  battery  as  a  linear  system  with  two  states,  namely,  Tco re  and  Tair, 
and  reference  [14]  uses  the  nonlinear  heat  transfer  equation  with  a 
single  state.  Simulations  show  that  the  dynamics  of  Tair  have 
negligible  fluctuations  around  the  ambient  temperature.  Therefore, 
the  temperature  model,  considered  below,  can  be  simplified  to  the 
linear  part  of  the  heat  transfer  equation 

T[k  +  1  ]  =  T[k]  -  a(T[k]  -  Tamb)  +  bi 2  [k]  (41 ) 


>'*['<]  k—  0, 1, ...,  It!  —  1  (37) 

It  is  seen  that  the  optimal  current  is  constant  and  is  a  function  of 
the  weight  on  TTC  and  the  series  resistance.  Therefore,  the  solution 
of  optimal  time-to-charge  and  energy  loss  (OtE)  problem  is  a 
CC— CV  profile  with  the  current  of  the  CC  stage  given  by  (37). 
Following  the  CC  stage,  from  Jq  to  kf,  one  has  the  CV  stage  where 


where 

A 

mbattOi,batt^Eff 

is  the  cooling  coefficient  and 


b  = 


RqA 

mbattQi.batt 


(42) 


(43) 


v[k\=vc  k  =  k1,...,kf  (38) 

To  the  best  of  our  knowledge,  this  is  the  first  time  that  it  is 
proved  that  the  well-known  CC— CV  charging  profile  is  the  optimal 
solution  of  a  particular  optimization  problem,  namely,  the  problem 
of  minimizing  the  weighted  sum  of  time-to-charge  and  energy  loss. 

In  the  sequel,  this  profile  is  referred  to  as  OtE  profile  or  OtE 
policy. 

Before  we  close  this  section,  we  point  out  another  way  of  solving 
the  OtE  problem  of  ( 12 )  and  ( 13 )  by  condensing  ( 13 )  for  all  values  of 
k  into  a  single  condition.  From  (13)  we  can  write 

i[k]  =  (s[k+\]-s[k])/ch  (39) 

Since  (39)  holds  for  k  =  0, 1, ....  /q  -  1,  and  using  the  initial  and 
end  values  of  SOC  from  (13)  we  can  write: 

ki  l 

J2  I'M  =  (Si  -  s0)/ch  (40) 

1=0 


Here  T  is  the  battery  core  temperature  in  kelvin  (K),  Tamb  is  the 
ambient  temperature  in  K,  mbatt  is  the  battery  mass  in  kg,  Ch  batt  is 
the  heat  capacity  of  the  battery  in  J/(kg  K),  and  REff  is  the  effective 
thermal  resistance  in  K/W  (kelvin/watt). 

Defining  temperature  rise  (TR)  as  T[k]  =  T[k]  -  Tamb  and 
assuming  T[0]  =  Tamb,  we  can  write 

t[k  +  1]  =  (1  -  a)t[k ]  +  bi2{k],  t [0]  =  0  (44) 

The  solution  of  (44)  is 

T{k]  =  bYl  0  —  o)k  1  li2[l]  (45) 

/=o 

Equation  (45)  states  that  the  temperature  rise  at  any  time  is  the 
integral  of  the  square  of  current,  from  time  zero  up  to  that  time 
with  a  “forgetting  factor”  of  (1  -  a)  and  the  scaling  factor  b. 

Since  T\k]  is  positive  for  any  k,  the  cost  function  including  TTC,  EL 
and  TR  can  be  written  as 


Therefore,  the  problem  of  (12)  and  (13)  is  equivalent  to  a 
quadratic  programming  problem  with  the  constraint  in  (40).  In  this 
way,  we  are  dealing  with  currents  i[/[  as  our  unknowns.  It  is  easy  to 
show  that  this  results  in  the  same  solution  as  (37).  This 


JtET  =  pjt  +  Je  +  PtJt  (46) 

where  Jt  and  JE  are  TTC  and  EL  as  before  and  JT  is  the  temperature 
rise  index  (TRI)  defined  as  follows: 
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;t  =  aE 

k= 0 

Since  f  [0]  =  0,  the  TRI  can  be  written  as 


(47) 


L  with  respect  to  A,  and  using  (56)  we  find  the  optimal  current 
profile  in  the  first  stage  as  follows: 


m  = 


geqW(£l  -Sp) 


k  =  o,i . It!  - 1 


(57) 


kf- 1 

h  =  T[k+  1]  (48) 

k= 0 


Using  (45)  and  (48),  we  can  write  (46)  as  follows 

kf  1 

JtET  =  Ptk/A  +  E  Ro!2[/<]A 
k—0 

kf- 1  k 

+PlbAj2  X)(1  —  a)fc-'i2[/]A 

k=0  1=0 


(49) 


which  can  be  simplified  as  follows 

%  i  /  %  k  '  \ 

JtET  =  Ptk/A  +  A  E  U>  +  Pt*>  E  (1  -  a)'  )  I2  [k]  (50) 

/<=o  \  i=o  ) 

Simplifying  the  inner  summation  and  noting  that  b/a  =  R0REff, 
we  can  write 

kf- 1 

JtET  =  ptkf  A  +  A  E  Req [k] I2  [k]  (51) 

k=0 


We  refer  to  the  current  profile  in  (57)  as  the  optimal  time-to- 
charge,  energy  losses  and  temperature  rise  (OtET)  policy.  Note 
that  (57)  is  similar  to  what  we  obtained  for  the  OtE  case.  In 
particular,  if  pT  =  0,  then  (57)  will  be  the  same  as  (28).  Also, 
comparing  (56)  with  the  OtE  case  and  noting  that  for 

k  =  0, 1 . ki  - 1  we  can  use  the  approximation  of 

Req[fc]=Ro(l  +PT^Eff)'  analogous  to  the  optimal  current  profile  of 
(37),  we  can  write 


m- 


p  t 


R0(l  +  PT^Eff) 


k  =  0, 1, ..., k1  —  1 


(58) 


We  refer  to  the  current  profile  of  (58)  as  the  near-optimal  time- 
to-charge,  energy  loss  and  temperature  rise  (NOtET)  policy. 


4.  Simulations 

In  this  section,  we  present  simulations  based  on  the  theoretical 
foundations  of  the  previous  section. 


4.3.  Verification  of  the  optimal  solution 


Req[k]  =R0  +  Ri[k]  (52) 

RT[k]  =  PTR0REff  (l  -  (1  -  a)kf-k )  (53) 

where  RT[k]  is  the  heating  equivalent  resistance.  Assume,  as  before, 
that  at  time  kf,  the  terminal  voltage  v  reaches  its  maximum 
allowable  value  of  vc,  and  SOC  reaches  S] .  Given  Sj  and  kf ,  we  can 
write  the  cost  function  as 


fci-i 

J(Sf,kf)  =  aE  Req[k]i2[k]  (54) 

k= 0 

Note  that  we  discarded  the  contributions  of  i[kf],  ■  ■■,i[kf  -  1], 
because  when  the  terminal  voltage  reaches  vc  the  current  is  already 
determined  by  the  constrained  dynamics  of  the  system  in  (9);  we 
also  discarded  the  contribution  of  kf,  i.e.  ptk/A,  because:  firstly,  kf  is 
given;  secondly,  given  Si,  kf  -  kf  is  also  known,  which  means  kf  is 
known.  An  important  point  to  note  is  that,  while  the  upper  bound 
of  the  summation  in  (54)  is  kf  -  1 ,  the  formulation  for  Req,  i.e.,  (52), 
considers  the  effect  of  the  whole  charging  time  and  it  contains  kf 
rather  than  kf . 

Now,  given  St  and  kf ,  we  can  state  the  optimal  charging  problem 
as  follows:  Minimize  (54)  subject  to  (40),  or  equivalently 


(k  i-i 

Minimize  :  L  —  J(Sf ,  kf )  +  A  E 

\t= o 


si  ^  s0\ 
ch  ) 


(55) 


Taking  the  derivative  of  Lagrangian  L  with  respect  to  i[k]  for  k  = 
0, 1, ...,  kf  —  1  and  equating  it  to  zero,  we  have: 


AGeq  [k] 
2A 


(56) 


2Req  [fc]  A 

where  Geq[k]  =  1  /Req[k]  is  the  conductance.  Taking  the  derivative  of 


Here,  we  apply  different  levels  of  current  and  the  simulation  is 
run  until  the  terminal  voltage  reaches  vc  and  after  that  a  constant 
voltage  of  vc  is  applied  until  the  battery  is  charged  to  skf.  Five 
different  current  profiles  are  chosen  including  the  optimal  current 
profile  (Fig.  3).  The  optimal  current  profile  as  mentioned  before  has 
the  value  of  y/pt/Ro  in  the  CC  stage.  The  battery  parameters  of 
Nokia  BP-4L  (Cell#3),  given  in  the  appendix,  are  used.  The 
following  simulation  parameters  are  used:  p  =  1,  A  =  l(s),  s0  =  0, 

Skf  =  1- 

The  appendix  also  shows  the  parameters  of  the  OCV  curve 
(calculated  based  on  [3]).  The  OCV  is  a  function  of  SOC  s  as  in  Ref. 

[3]. 

z(s)=E  +  s(l  -  2E)  (59) 


OCV  (z)  =  K0  +  KfZ -1  +  I<2z-2  +  I<3z-3  +  I<4z-4  +  I<5z  +  K6\n{z) 


+  /C7ln(l  -z) 


(60) 


and  £  =  0.15.  Fig.  3  shows  the  current  profiles  with  different 
levels  of  current  in  the  CC  stage.  As  seen  from  Fig.  3,  at  lower  levels 
of  current,  the  CC  stage  will  take  a  longer  time  and  the  terminal 
voltage  reaches  the  threshold  voltage  of  vc  at  a  later  time.  At  higher 
levels  of  current,  however,  the  OCV  grows  more  rapidly.  As  the 
terminal  voltage  is  v[k]  =  V0  (s[k])  +  R0i[k],  at  higher  levels  of  cur¬ 
rent  the  threshold  voltage  of  vc  is  reached  in  a  shorter  time. 

Fig.  4(a)  shows  the  cost  function  JtE  for  the  five  current  profiles 
of  Fig.  3  and  Fig.  4(b)  shows  the  corresponding  current  levels  in  the 
CC  stage.  It  is  seen  that  the  optimal  current  profile  (i.e.,  profile  3) 
has  the  lowest  cost  function.  Deviating  from  this  profile,  either  by 
increasing  or  decreasing  the  current  in  the  CC  stage,  results  in  an 
increase  in  the  cost  function.  For  the  lower  current  levels  (profiles 
1—2),  the  rise  in  the  cost  function  is  due  to  a  rise  in  TTC  and  for 
higher  current  levels  (profiles  4—5)  the  rise  in  cost  function  is  due 
to  rise  in  EL. 
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Fig.  3.  Five  different  current  profiles  (including  the  optimal  profile). 


(a) 


r 


1  2  3  4  5 

profile  number 
(b) 


1  2  3  4  5 

profile  number 

Fig.  4.  Five  different  current  profiles:  a)  cost  functions  b)  current  levels  in  CC  stage. 


4.2.  Effect  of  weights 

In  this  subsection,  we  use  different  cost  functions  and  find  the 
corresponding  optimal  profiles.  Different  values  of  pt  from  0.1  to  0.5 
are  chosen.  Figs.  5—7,  respectively,  show  the  profiles  of  current, 
state  of  charge  and  terminal  voltage.  Fig.  5  shows  that  low  values  of 


pt  result  in  low  values  of  current  in  the  CC  stage.  In  other  words,  a 
low  pt  puts  less  emphasis  on  charging  time  and  more  emphasis  on 
the  energy  losses;  hence,  it  results  in  low  level  of  current  which 
provides  low  energy  losses.  On  the  other  hand,  by  increasing  pt, 
more  emphasis  is  placed  on  the  charging  time.  Consequently,  the 
level  of  current  is  increased  proportionally  to  ^fp[  to  reduce  the  TTC. 


Fig.  5.  Current  profiles  for  different  values  of  pt. 
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Fig.  6.  SOC  profiles  for  different  values  of  pt. 


Fig.  7.  Voltage  profiles  for  different  values  of  pt. 


Fig.  6  shows  the  state-of-charge  profiles  for  different  values  of 
pt.  It  is  seen  that  by  increasing  pt,  more  emphasis  is  placed  on 
charging  time  and  the  SOC  reaches  the  final  value  in  a  shorter  time. 
Fig.  7  shows  the  terminal  voltage  profiles  for  different  values  of  pt. 
Note  that  for  low  values  of  pt,  as  the  emphasis  on  energy  loss  is 
high,  the  corresponding  current  level  in  the  CC  stage  is  low,  and 
consequently,  the  terminal  voltage  reaches  the  threshold  value  of  vc 
at  a  later  time.  Hence,  the  duration  of  the  CC  stage  is  high  and  the 
charging  time  is  high  as  well. 

Fig.  8  shows  the  time-to-charge,  energy  losses  and  efficiency  as 


functions  of  pt.  As  expected,  high  values  of  pt  result  in  lower  TTC. 
The  low  TTC,  however,  is  obtained  by  increasing  the  current  level; 
as  EL  is  proportional  to  the  square  of  current,  thus  the  high  values  of 
pt  result  in  high  values  of  EL.  The  high  values  of  EL  mean  that  a 
higher  fraction  of  input  power  is  wasted;  hence  it  is  equivalent  to  a 
decline  in  efficiency. 

Fig.  9  shows  the  time-to-charge  versus  efficiency  (ratio  of 
effective  to  total  energy)  curve.  TTC  and  efficiency  are  two  coun¬ 
teracting  objectives.  For  low  values  of  pt,  as  less  emphasis  is  put  on 
TTC,  the  TTC  is  high;  however,  high  TTC  is  the  result  of  low  current 


Pt 


Fig.  8.  TTC,  EL  and  efficiency  curves  for  different  values  of  pt. 
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Table  1 

Battery  thermal  parameters. 


Paramter  set 

mbatt  (kg) 

REff  (K/W) 

Qi.batt  U/(kg-I<)) 

A 

0.37824 

7.8146 

795 

B 

0.080 

1.6528 

168.15 

Table  2 

Cost  function  for  different  schemes. 

Pt 

Pt 

Thermal  parameters 

JtET 

OtE 

NOtET 

OtET 

1 

1 

A 

26,734 

22,206 

22,198 

1 

4 

A 

72,023 

38,996 

38,960 

1 

1 

B 

14,970 

14,696 

14,696 

1 

4 

B 

24,966 

21,272 

21,271 

values,  which  incur  low  energy  losses  and  hence  higher  efficiency. 
For  example  at  pt  =  0.1,  the  TTC  is  195  min,  but  the  efficiency  is  as 
high  as  95.82%.  On  the  other  hand,  for  high  values  of  pt  which  place 
more  emphasis  on  TTC,  the  TTC  is  reduced  dramatically;  however, 
low  TTC  is  achieved  by  increasing  the  current  values,  which  results 
in  high  energy  losses  and  hence  lower  efficiency.  For  example,  at 
pt  =  0.5,  the  TTC  is  as  low  as  148  min,  but  the  efficiency  decreases  to 
92.87%. 

4.3.  Temperature  effect 

In  this  section,  we  consider  the  effect  of  temperature  rise  index 


(TRI)  on  optimal  charging.  The  cost  function  is  a  weighted  sum  of 
TTC  (seconds),  EL  (Joules)  and  TRI  (Kelvin  seconds),  given  by 

JtET  =  Pt  x  TTC  +  EL  +  pT  x  TRI  (61 ) 

We  used  two  sets  of  thermal  parameters,  shown  in  Table  1. 
Parameter  set  “A”  is  adopted  from  Ref.  [19],  Parameter  set  “B"  is  a 
scaled  version  of  parameter  set  “A”  with  mbatt  set  as  the  weight  of 
Nokia  BP-4L.  For  each  set  of  thermal  parameters  (“A”  or  “B”),  the 
weights  of  the  cost  function  are  chosen  as  pt  =  1,  px  =  1  and  pt  =  1, 
pT  =  4.  Three  schemes  are  used:  OtE  (equations  (37)  and  (38)),  OtET 
(equations  (38)  and  (57)),  and  NOtET  (equations  (38)  and  (58)).  The 
cost  function  in  (61)  or  (46)  is  calculated  for  the  three  schemes. 
Table  2  shows  the  cost  functions  of  the  three  schemes  for  different 
weightings.  As  seen  from  this  table,  the  cost  function  for  the  OtE 
has  the  highest  value.  Also  the  difference  between  the  cost  function 
of  OtET  and  NOtET  is  negligible  with  the  OtET  being  slightly  smaller 
when  thermal  parameter  set  “A”  is  used.  For  thermal  parameter  set 
of  “B”,  there  is  visually  no  difference  between  NOtET  and  OtET.  Due 
to  this  negligible  difference  in  the  cost  function  and  also  since  the 
calculation  of  NOtET  profile  is  much  easier  than  that  of  the  OtET,  it 
is  reasonable  to  use  NOtET  rather  than  the  OtET  scheme.  Also  note 
that  the  weight  on  TRI  results  in  a  reduction  of  current,  as  can  be 
seen  from  Fig.  10.  This  reduction  in  current  level  results  in  a  lower 
temperature  rise  (see  Fig.  11 ).  In  other  words,  energy  losses  with  Req 
instead  of  R0  can  be  used  as  a  surrogate  cost  function  for  the  TRI. 

4.4.  Analysis  of  different  commercial  batteries 

In  this  section,  we  discuss  the  behavior  of  different  commercial 
batteries.  The  parameters  of  the  investigated  batteries  are  given  in 
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time  (min) 

Fig.  11.  Temperature  profiles  for  pt  =  l,pT  =  1  and  temperature  parameter  set  “A”. 
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Fig.  12.  Time  to  charge  versus  efficiency  of  different  battery  types  at  25  °C. 


the  appendix.  Note  that  the  equivalent  electrical  circuit  parameters 
given  in  Table  3  are  for  model  111  (see  Fig.  2).  In  simulations,  we  use 
the  summation  of  “R0  +  R of  model  III  as  an  estimate  of  resistance 
R0  in  model  I.  The  batteries  are  Samsung  EB575152  (four  cells), 
Samsung  EB504465  (four  cells),  Samsung  AB463651  (two  cells), 
Nokia  BP-4L  (four  cells),  LG  LGIP  (two  cells). 

Next,  we  apply  the  OtE  algorithm  with  pt  =  0.5  to  16  commercial 
batteries  to  investigate  the  times-to-charge  and  efficiencies  of  the 
batteries.  The  parameters  of  the  batteries,  i.e.,  the  electrical 


parameters  of  the  models  in  Figs.  1  and  2,  and  the  parameters  of  the 
OCV  function  in  (60),  were  calculated  using  experimental  data  and 
by  applying  the  BFG  algorithms  in  Refs.  3,4],  These  parameters  are 
listed  in  Tables  3  and  4  in  the  Appendix.  Fig.  12  shows  the  TTC 
versus  efficiency  for  different  types  of  batteries.  Among  all  batte¬ 
ries,  Sam-EB575152  (Cell  3)  has  the  lowest  efficiency  (90.73%).  This 
can  be  attributed  to  the  high  resistance  of  this  battery,  which  might 
be  due  to  aging.  Sam-EB504465  (Cell  4)  has  the  highest  TTC 
(102  min)  and  Nokia  BP-4L  (Cell  4)  has  the  highest  efficiency.  Note 


Fig.  13.  Cost  function  for  different  battery  types  at  25  °C,  pt  =  0.5. 
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that  the  cells  of  the  same  battery  are  close  to  each  other  in  terms  of 
efficiency  and  TTC.  Considering  all  the  cells  of  a  battery,  we  can  say 
that  LG-LGIP  cells  (circle  markers)  have  the  highest  efficiency 
(91.4%).  Fig.  13  shows  the  cost  function  values  of  JtE  =  Pt/t  +  h- 
When  TTC  is  weighted  with  weight  value  of  pt  =  0.5,  Sam- 
EB575152  (Cell  2)  has  the  best  performance. 

5.  Conclusion 

The  optimal  charging  problem  involving  a  weighted  combina¬ 
tion  of  time-to-charge  (TTC),  energy  loss  (EL)  and  temperature  rise 
index  (TRI)  was  considered.  The  optimal  TTC  and  EL  solution  (OtE) 
is  found  to  be  the  well-known  CC— CV  strategy  with  the  value  of 
current  in  the  CC  stage  being  a  function  of  the  ratio  of  weighting  on 
TTC  and  EL  and  also  the  resistance  of  battery.  To  the  best  of  our 
knowledge,  this  is  the  first  time  that  it  is  proved  that  the  well- 
known  CC— CV  charging  profile  is  the  optimal  solution  of  a  partic¬ 
ular  optimization  problem,  namely,  the  problem  of  minimizing  the 
weighted  sum  of  time-to-charge  and  energy  loss.  In  addition,  an 
analytical  solution  for  the  optimal  TTC,  EL  and  TRI,  referred  to  as 
OtET,  was  developed.  Due  to  similarity  of  the  structure  of  the  OtE 
and  OtET  solutions,  a  near-optimal  version  of  OtET  was  developed 
(referred  to  as  NOtET).  The  NOtET  is  a  CC— CV  strategy  with  the 
value  of  current  in  the  CC  stage  being  a  function  of  the  ratio  of 
weighting  on  TTC  and  EL,  the  resistance  of  the  battery  and  the 
effective  thermal  resistance.  A  number  of  simulations  were  con¬ 
ducted  to  evaluate  the  effect  of  weighting  parameters.  Finally, 
extensive  results  on  industrial  batteries  from  LG,  Nokia  and  Sam¬ 
sung  were  presented. 
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Nomenclature 

a,  b  coefficients  in  temperature  model,  (42,43) 

Cbatt  capacity,  (2) 

ch  parameter  in  coulomb  counting,  (1) 

Qj  batt  heat  capacity  of  the  battery  in _//(kg-K),  (42) 

CC  constant  current 

CV  constant  voltage 

A  sampling  time,  (2) 

EL  energy  loss 

t 7  charging  efficiency 

Geq[k]  conductance  equal  to  reciprocal  of  Req[k],  (56) 

H[k]  Hamiltonian  function,  (16) 

i[k]  charging  current,  (1) 

t*[k]  optimal  current,  (28) 

i*  [fc|/<i  ]  optimal  current  given  /<q 

Je  energy  loss  cost  function,  (3) 

Jj  temperature  rise  cost  function,  (46) 

Jt  time  to  charge  cost  function,  (3) 

JtE  objective  function  as  a  combination  of  TTC  and  EL  cost 

functions,  (12) 

/tE*  partially  optimized  cost  function,  (30) 

/tEX  cost  function  including  time-to-charge  (TTC),  energy  loss 
(EL),  and  temperature  rise  index  (TRI),  (46) 


JE(ki)  minimum  energy  loss  cost  function  as  a  function  of  kq 
(given  ki ),  (29) 

/*E(kj)  optimized  cost  function 

k  time  index,  (1) 

k\  optimal  final  time  for  constant-current  charging  stage, 

(34) 

k,+  solutions  to  partially  optimized  cost  function,  (33) 

/q  time  when  terminal  voltage  reaches  maximum  allowed 

value,  (11) 

kf  final  value  of  time  index,  (3) 

L  Lagrangian  function,  (55) 

A[k]  co-state  variable,  (18) 

mbatt  battery  mass  in  kg,  (42) 

NOtET  near  optimal  TTC,  energy  loss,  and  temperature 

OCV(z)  open  circuit  voltage,  (60) 

OtE  optimal  time-to-charge  and  energy  loss 
OtET  optimal  TTC,  energy  loss,  and  temperature  rise 
R0  resistance  in  model  I  of  battery  as  shown  in  Fig.  1,  (3) 

REff  effective  thermal  resistance  in  K/W,  (53) 

Req[k]  heating  equivalent  resistance  added  to  the  battery 
resistance,  (52) 

Rx[k]  heating  equivalent  resistance,  (52) 

pt  ratio  of  weight  of  TTC  cost  function  to  weight  of  EL  cost 

function,  (12) 

s  state  of  charge,  (1) 

s0  initial  SOC,  (4) 

Sj  state  of  charge  when  terminal  voltage  reaches  maximum 

allowed  value,  (11) 
skf  final  SOC,  (4) 

T  battery  core  temperature  in  kelvin,  (41 ) 

((mil)  ambient  temperature  in  kelvin,  (41 ) 

T[k]  temperature  rise,  (44) 

TTC  time-to-charge 

v[k]  terminal  voltage  of  battery,  (38) 

V0  open-circuit  voltage,  (5) 

vc  voltage  corresponding  to  SOC  of  1,  (8) 

wE  weight  for  energy  loss  cost  function,  (3) 

wt  weight  for  time  to  charge  cost  function,  (3) 

z(s)  parameter  relating  OCV  to  SOC,  (59) 


Appendix 

The  following  tables  show  the  parameters  of  the  equivalent 
electrical  circuit  model  HI  for  different  commercial  batteries. 


Table  3 

Electrical  parameters  of  model  III  for  commercial  batteries. 


Make 

Model 

Cell# 

R0  (mfi) 

R !  (mC2) 

Q  (F) 

a 

Cbatt  (Ah) 

Samsung 

EB575152 

1 

253 

106 

4581 

0.997934 

1.1875 

Samsung 

EB575152 

2 

209 

94 

5203 

0.997962 

1.2187 

Samsung 

EB575152 

3 

418 

58 

6222 

0.99724 

1.2001 

Samsung 

EB575152 

4 

200 

142 

3046 

0.997689 

1.485 

Samsung 

EB504465 

1 

259 

106 

4598 

0.997941 

1.5001 

Samsung 

EB504465 

2 

268 

168 

2493 

0.997615 

1.5293 

Samsung 

EB504465 

3 

272 

211 

1680 

0.997186 

1.5261 

Samsung 

EB504465 

4 

287 

224 

1589 

0.997189 

1.4831 

Samsung 

AB463651 

1 

451 

198 

2100 

0.997597 

0.9791 

Samsung 

AB463651 

2 

294 

214 

1950 

0.997602 

0.9614 

Nokia 

BP-4L 

1 

263 

100 

5031 

0.998012 

1.5514 

Nokia 

BP-4L 

2 

264 

64 

8141 

0.99808 

1.5691 

Nokia 

BP-4L 

3 

258 

95 

5306 

0.998028 

1.5612 

Nokia 

BP-4L 

4 

228 

50 

10502 

0.998106 

1.613 

LG 

LGIP 

1 

264 

101 

4747 

0.997919 

1.1141 

LG 

LGIP 

2 

297 

76 

6654 

0.998021 

1.1121 
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Table  4 

OCV  parameters  for  commercial  batteries. 


Make 

Model 

Cell# 

OCV  parameters 

Ko 

*1 

k2 

'<3 

«4 

k5 

K6 

Kj 

Samsung 

EB575152 

1 

-1.1288 

65.0931 

-10.5332 

1.0640 

-0.0457 

-52.2654 

94.1916 

-0.7417 

Samsung 

EB575152 

2 

-1.6976 

68.3306 

-11.1189 

1.1290 

-0.0487 

-54.5366 

98.4237 

-0.7882 

Samsung 

EB575152 

3 

-1.0802 

68.3600 

-11.0617 

1.1178 

-0.0480 

-55.3010 

99.0817 

-0.8059 

Samsung 

EB575152 

4 

-4.3113 

22.9007 

-4.2921 

0.4926 

-0.0239 

-10.3460 

27.2251 

0.0226 

Samsung 

EB504465 

1 

0.0218 

54.9000 

-9.1299 

0.9668 

-0.0442 

-44.6447 

78.9622 

-0.7934 

Samsung 

EB504465 

2 

1.8254 

61.0951 

-10.0031 

1.0470 

-0.0474 

-52.6013 

89.7801 

-1.0640 

Samsung 

EB504465 

3 

0.2510 

55.0370 

-9.1256 

0.9634 

-0.0439 

-45.0942 

79.3945 

-0.8297 

Samsung 

EB504465 

4 

2.9648 

59.8808 

-9.7283 

1.0127 

-0.0457 

-52.9562 

88.9546 

-1.1311 

Samsung 

AB463651 

1 

-1.6972 

41.8528 

-7.0700 

0.7522 

-0.0343 

-30.6508 

58.1983 

-0.4098 

Samsung 

AB463651 

2 

-1.2526 

40.3216 

-6.7711 

0.7166 

-0.0326 

-29.7536 

56.3932 

-0.3814 

Nokia 

BP-4L 

1 

-3.2203 

51.9246 

-8.8187 

0.9344 

-0.0421 

-38.1050 

71.7162 

-0.5991 

Nokia 

BP-4L 

2 

-2.7537 

52.9707 

-8.9327 

0.9407 

-0.0422 

-39.6357 

73.7620 

-0.6418 

Nokia 

BP-4L 

3 

-3.2084 

51.8554 

-8.7993 

0.9314 

-0.0419 

-38.0572 

71.6483 

-0.5996 

Nokia 

BP-4L 

4 

-2.7140 

60.3626 

-10.0810 

1.0533 

-0.0469 

-46.2542 

85.0092 

-0.7139 

LG 

LGIP 

1 

0.5267 

61.5448 

-10.1553 

1.0682 

-0.0485 

-51.2165 

89.3849 

-0.9091 

LG 

LGIP 

2 

0.4788 

59.0975 

-9.7677 

1.0290 

-0.0468 

-48.9737 

85.6643 

-0.8748 
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